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ABSTRACT 

This  paper  presents  a  computational  investigation  of  the  validity  of  eddy  current  testing  (ECT)  for  defects  embedded  in  steel  using 
parametrically  designed  defects.  Of  particular  focus  is  the  depths  at  which  defects  can  be  detected  through  ECT.  Building  on  this  we 
characterize  interior  defects  by  parametrically  describing  them  and  then  examining  the  response  fields  through  measurement.  Thereby 
we  seek  to  establish  the  depth  and  direction  of  detectable  cracks.  As  a  second  step,  we  match  measurements  from  eddy  current 
excitations  to  computed  fields  through  finite  element  optimization.  This  develops  further  our  previously  presented  methods  of  defect 
characterization.  Here  rough  contours  of  synthesized  shapes  are  avoided  by  a  novel  scheme  of  averaging  neighbor  heights  rather  than 
using  complex  Bezier  curves,  constraints  and  such  like.  This  avoids  the  jagged  shapes  corresponding  to  mathematically  con  ect  but 
unrealistic  synthesized  shapes  in  design  and  nondestructive  evaluation. 

CITATION:  Mathialakan.  T..  Kartliik.  V.,  Jayakumar.  P..  Thyagarajan.  R.  et  al..  "The  Depth  Limits  of  Eddy  Current  Testing  for  Defects:  A 
Computational  Investigation  and  Smooth-Shaped  Defect  Synthesis  from  Finite  Element  Optimization."  SAEInt.  J.  Mater.  Manf.  8(2):2015. 
doi:l 0.427 1/20 15 -01 -05 95. 


INTRODUCTION 

Eddy  current  testing  (ECT)  is  a  widely  accepted,  cheap  and  portable 
method  for  the  detection  of  cracks  and  other  defects  in  conductive 
materials.  While  ECT  is  a  favorite  with  engineers  for  surface  defect 
identification  through  nondestructive  evaluation  (NDE).  for  more 
deeply  embedded  defects  the  literature  is  not  clear.  It  lias  been 
established  that  low  frequency  testing  can  distinguish  “all  the 
different  levels"  of  corrosion  w  hen  a  sample  is  exposed  to  a  very 
corrosive  environment  for  3  months  QJ  while  other  researchers  have 
established  that  eddy  current  testing  can  identify  corrosion  to  depths 
down  to  15  mm  [2.  3J.  Submarines  and  aircraft  employ  eddy  current 
probes  for  testing  thin  (up  to  1  inch)  conducting  systems  [4j.  and 
ultrasound  [5]  for  thicker  submarine  hulls  [6J.  For  deeper  penetration, 
pulsed  eddy  currents  [7]  or  remote  field  eddy  currents  are  used  [8.  9]. 
ECT  Technology  is  therefore  perfect  for  the  piupose  at  hand. 


Flg.l.  Penetration  of  Eddy  Currents  in  the  Presence  of  Defects 


However,  the  reality  is  that,  thinking  qualitatively  about  the  three 
situations  in  Fig.l.  the  depth  for  detection  depends  on  defect 
orientation.  In  the  first,  when  a  defect  is  parallel  to  the  surface,  the 
defect  will  be  detectable  through  the  interruption  of  the  eddies.  In 
the  second  situation,  it  is  perpendicular  to  the  surface  but  not 
reaching  the  surface:  we  would  expect  the  skin-effect  to  allow  the 
currents  to  flowT  along  the  surface.  So  it  will  be  hard  to  detect  the 
defect  in  the  second  case.  Whereas  in  the  third  situation  if  the 
defect  is  perpendicular  to  the  surface  w  ithout  a  gap  for  the  eddy 
current  to  flowr  along  the  surface,  the  defect  w'ould  interrupt  the 
flowr  of  eddy  currents  w  hich  wrould  find  it  difficult  to  go  deep 
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because  of  the  skin  effect.  Therefore  this  will  be  easily  detectable. 
On  the  other  hand  if  the  defect  is  thin  and  long,  a  3D  model  will 
show  no  current  interruption. 


where  //  is  the  spatial  permeability  distribution  and  a  is  the 
conductivity.  The  magnetic  flux  density  B  is  related  to  the  magnetic 
vector  potential  as, 


This  interruption  would  make  the  presence  of  defects  easily 
detectable.  Different  types  of  defect  models  with  the  ECT  coil  are 
shown  in  Fig.2.  Defects  a,  b,  c,  d  &  e  in  Fig.2  show  the  difference  in 
angle,  depth  and  length.  Moreover,  if  the  axis  of  the  external  exciting 
ECT  coil  is  parallel  to  the  surface,  our  considerations  need  to  be 
entirely  different.  This  study  therefore  through  a  series  of  finite 
element  analyses  of  parameterized  geometries  of  a  single  defect, 
seeks  to  identify  the  limits  of  eddy  current  testing  in  NDE.  This  study 
is  confined  to  steel  plates  for  army  ground  vehicle  armor.  A 
qualitative  relationship  between  the  depth  of  a  detectable  defect,  the 
frequency  of  excitation  and  the  geometric  parameters  is  sought  to 
help  engineers  choose  whether  or  not  to  use  ECT. 


Fig.2.  Defect  Model:  The  model  has  defects  in  different  angle  (a,  b,  c),  depth 
(c,  e),  and  length  (d,  e). 


DETECTION  OF  DEFECT:  FIELD  CHANGES 

A  defect  in  a  material  can  be  varied  in  shape  and  location.  The  defect 
detection  method  should  concentrate  on  the  characteristics  of  the 
defect.  These  effects  are  analyzed  in  this  paper  for  varying  size,  shape 
and  location  of  defect  in  materials.  An  iterative  approach  is  presented 
that  repeatedly  employs  the  finite  element  technique  for  modeling  the 
forward  problem  to  calculate  the  effect  caused  by  the  defects  in  a 
steel  plate. 


V  x  A  =  B 


(2) 


In  two  dimensional  problems,  the  current  density  J  and  the  magnetic 
vector  potential  A  will  be  in  the  z  direction  (i.e.  the  transverse  direction 
in  which  no  changes  occur).  Therefore  the  magnetic  field  density  B  will 
be  in  the  X  and  /  directions.  From  (1)  and  (2)  we  can  write, 


±1 

,  ±1 

f— 'I 

ax 

\[ldx  J 

dy 

\fidyJ 

(3) 


Finite  element  analysis  [10]  provides  the  solution  to  (1)  by 
applying  certain  boundary  conditions.  This  leads  to  the  finite 
element  matrix  equation 


[P]MJ  =  {R} 


(4) 


where  [P]  is  the  finite  element  stiffness  matrix  and  {R}  is  a 
column  vector. 


Any  defect  in  the  material  should  affect  the  magnetic  field  density  B 
between  the  AC  coil  and  the  material.  The  magnetic  field  B  is 
calculated  at  a  line  called  the  measuring  line  which  is  divided  into  m 
points.  The  effect  in  B  caused  by  the  defect  for  the  m  measuring 
points  can  be  calculated  using  the  expression 

D0  _  irinY  ( \Bnodefect ~  gde/ect|\ 

/W)  -  maxi=0..m  (J  BU/ec[  |  j 

(5) 

rB 

(O.d.i)  is  the  maximum  value  of  flux  ratio  between  flux  changes  and 
flux  without  the  defect,  where  a  defect  with  length  l  is  at  depth  d  from 
the  material  surface  and  rotated  by  angle  6  clock- wise  from  the 
“horizontal”  line  parallel  to  the  steel  surface. 


We  calculate  the  magnetic  field  density  for  the  known  defects  from 
the  forward  problem  using  the  finite  element  method  [2]  and  examine 
the  extent  to  which  the  exterior  field  is  altered  by  the  defect  to  see  if 
the  presence  of  the  defect  can  be  discerned  through  the  changes  in  the 
measurable  external  field. 

PROBLEM  STATEMENT 

Magnetic  fields  in  a  ferromagnetic  material  can  be  generated  by 
placing  an  AC  (Alternative  Current)  coil  on  top  of  the  material.  For 
AC  magnetization,  the  vector  magnetic  potential  A,  and  the  exciting 
current  density  J  at  angular  frequency  co,  are  related  by 

-V  x  -  V  x  A  =  I  -  jrnoA 

u 


METHODOLOGY 


Fig.3.  Design  Cycle  for  the  computation  process 


(1) 
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The  computational  process  in  the  defect  identification  system  is 
shown  in  Fig. 3.  It  calculates  the  effect  on  the  flux  density  for 
different  kinds  of  defect  that  vary  in  their  characteristics:  depth  (< d ), 
angle  ( 6 )  and  length  (/).  The  characteristics  of  a  defect  are  the 
changing  parameters  in  each  calculation  of  the  above  computational 
process.  In  the  latter  part  of  our  work,  we  extend  this  model  for 
design  optimization.  Mesh  generation  is  a  very  important  part  of 
finite  element  analysis  based  design  optimization.  We  need  to  use  a 
parameter  based  mesh  generator,  because  in  each  iteration  the  mesh 
has  to  be  generated  automatically  when  design  parameters  change. 
Therefore  we  use  a  special  script-based  parametric  mesh  generator 
HU  using  as  backend  the  single  problem  mesh  generator  Triangle 
[12]  to  get  the  corresponding  finite  element  solution  for  the 
magnetic  vector  potential  A.  Then  from  A,  we  compute  the 
magnetic  field  density  B.  From  B  we  calculate  the  effect  on  the  flux 

rB 

density  In  each  iteration,  the  electric  field  map  is  also 

generated  and  analyzed  visually. 

APPLICATION 

The  numerical  example  in  Fig.4  is  used  to  validate  the  proposed 
algorithm.  The  coil  (with  n=  1.0,  current  density  J  =  ±  1  A/m2,  and  co 
=  10  rad/s)  excites  the  magnetic  field  in  the  steel  plate  (with  jur  = 
100.0  and  current  density  J  =  0.0).  The  conductor  is  surrounded  by  air 
(with  jur=  1.0  and  current  density  J=  0.0).  The  magnetic  field  density 
in  the  V  direction  By  is  calculated  at  1 0  points  as  shown  in  Fig.4. 
labeled  as  the  measuring  line. 

First,  By  for  the  steel  plate  with  no  defect  is  calculated  at  the 
measuring  line  and  named  Bnodefect.  Since  our  study  is  about  whether 
the  defect  is  detectable  or  not,  we  defined  a  single  defect  with 
parameterized  geometries  {x,  d,  /,  w,  0}  (Fig.4.). 


Fig.4.  Numerical  Model 

By  changing  the  parameters  of  the  defect,  we  change  the  position  of 
the  defect  as  shown  in  Fig.5.  and  calculate  By  at  the  measuring 
points,  which  we  name  Bdefect.  Then  the  effect  on  flux  density 
R(d,d,i)  was  calculated  for  each  defect  and  tabulated.  By  considering 
the  maximum  of  R,  we  establish  the  defect  detecting  quantity  which 
is  measured. 


Fig.5.  Equipotential  lines  for  the  defects  at  various  positions 


Mathialakan  et  al  /  SAE  Int.  J.  Mater.  Manf.  /  Volume  8,  Issue  2  (May  2015) 


d.  d  =  0.3  and  0  =  90 


Fig.5.  (cont.)  Equipotential  lines  for  the  defects  at  various  positions 


Ratio  with  depth  for  horizontal  defect 


a. 


W,0 


varies  with  d  for  horizontal  defect 


Fig-6. 


Ratio  with  angle 


rb 

c  (d.d.l)  varies  with  0 

Rb 

Fig.6.  (cont.)  with  defect  location 

rB 

Fig  6  shows  how  the  ratio  varies  with  defect  locations.  When 

(O.d.i)  is  high,  there  is  a  higher  chance  that  the  defect  will  be 
detected.  Fig  6. a  shows  the  results  for  a  horizontal  defect  (0  =  0°.) 
when  the  depth  d  increases:  (0«d»O  decreases.  So  that  means  when 
the  depth  of  the  defect  is  higher,  it  would  be  harder  to  detect  as  to  be 
pvrwted.  Fig  6.b  shows  when  the  defect  is  vertical  (0  =  9°  ^  how 
(0.d.O  varies  with  depth.  Here  also  when  d  increases  (0-d«O 
decreases.  Here  d  =  0  corresponds  to  the  right-most  part  of  Fig.l. 
Anyway  when  we  compare  a  horizontal  defect  and  vertical  defect,  the 
horizontall  defect  has  a  higher  chance  to  be  detected.  This  is  to  be 
expected  because  a  horizontal  defect  will  interrupt  fluv^nnre  since 

the  flux  flow  i«  glnncr  the  surface.  Fig.6. c  shows  how  varies 

R  ^ 

with  angle  0:  (0-d*O  goes  to  a  maximum  when  0  =  45°.  So  we  could 

say  that  when  the  angle  of  the  defect  is  45°,  there  is  higher  chance  to 
be  detected. 

Experimental  work  was  done  to  verify  the  computational  work  of 
defect  detection.  The  coil  was  moved  along  the  x  axis  of  the  steel 
plate  (Fig.7.b.)  and  the  voltage  induced  was  measured  as  shown  in 
Fig/7.  We  created  a  rectangular  defect  on  a  steel  plate.  As  shown  in 
Fig.7.a.  A  close-up  look  of  the  defect  is  shown  in  Fig.7.b.  Alternative 
sinusoidal  current  was  excited  in  the  coil  with  300  mV  at  20  kHz 
using  the  waveform  generator.  A  pick  up  coil  was  wound  on  top  of 
the  AC  coil  to  measure  the  voltage  induced.  Since  the  voltage 
induced  is  very  small,  we  used  the  lock-in  amplifier  to  measure 
voltage.  We  moved  the  coil  along  the  surface  of  the  steel  and 
measured  the  induced  voltage  at  the  measuring  points  indicated  on 
Fig/7.  According  to  Faraday’s  equation,  the  voltage  induced  Vind  is 
proportional  to  the  change  in  flux  (p  (6). 

The  results  were  tabulated  and  plotted  in  Fig.  8.  From  the  plot,  when 
the  coil  is  just  above  the  defect,  the  voltage  induced  is  at  its 
maximum.  The  test  was  repeated  with  different  input  voltages  and 
different  frequency  and  we  obtainedt  the  same  behavior  as  Fig.8.  In 
this  way  we  could  detect  the  defects  and  moved  to  the  defect 
characterization  part  of  this  study.  The  frequency  of  excitation  also 
plays  an  important  role  in  detecting  the  defect.  As  in  equation  (6), 
when  frequency  co  goes  high,  Vind  will  be  high  and  can  be  measured 
easily.  But  the  eddy  current  will  flow  close  to  the  surface  due  to  the 
skin  effect.  Therefore  we  cannot  detect  defects  which  are  deeply 
inside  the  plate.  Our  lock-in  amplifier  (Model  SR844)  has  a 
bandwidth  of  20  kHz  to  200  MHz.  Due  to  the  limitation  of 
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resources;  we  could  not  experiment  with  low  frequencies.  We  did  yind 

our  work  using  20  kHz,  the  lowest  frequency  that  we  can  use  to 
measure  the  readings. 


EXPERIMENTAL  WORK 


Fig. 7.  Experimental  setup 

ECT  Probe 


Waveform 

Generator 


Lock-In 

Amplifier 


Steel  Plate 


Fig.7.a.  Lab  setup 


Fig.7.b.  Close-up-look  of  the  defect 


DEFECT  CHARACTERIZATION 

After  detecting  the  defect,  we  investigated  more  on  defect 
characterization.  It  is  important  to  know  the  size  or  character  of 
the  defect  after  establishing  that  there  is  defect  inside.  So  in  our 
work  we  investigate  and  establish  a  procedure  for  defect 
characterization  so  that  a  decision  to  withdraw  a  defective  part 
may  be  thought-out,  justifiable. 


AC 

o  o 


The  methodology  at  present  [13]  examines  the  response  of  the  hull 
under  test  to  an  excitatory  signal  from  an  eddy  current  probe.  By 
knowing  the  response  when  there  is  no  defect,  if  the  response  is 
different  because  of  the  defect,  the  test  object  is  presently  flagged  as 
defective  and  the  plate  is  sent  for  repairs  without  assessing  if  the 
defect  is  serious  enough  for  removal  from  service.  In  our  work,  we 
extend  that  methodology  to  defect  characterization.  An  iterative 
approach  is  presented  that  repeatedly  employs  the  finite  element 
technique  for  modeling  the  forward  problem  to  characterize  the  shape 
of  defects  in  a  steel  plate. 

We  can  calculate  the  magnetic  field  density  for  the  known  defects 
from  the  forward  problem.  But  in  our  inverse  problem  we  need  to 
know  the  characteristics  of  the  defect  for  that  field  configuration.  In 
design  optimization,  the  problem  geometry  is  defined  in  terms  of 
design  parameters  contained  in  a  vector  ^  =  (  x\»  3^1  >  yni 
(Fig.9.).  An  objective  function  F  is  defined  as  the  sum  of  the  squares 
of  the  difference  between  computed  and  measured  (defect) 
performance  values:  at  measurement  points  i, 

y\>  •  •  t y?i)  ~  2i(®Calculated  ~  ^Measured  ) 

(5) 

F  is  a  function  of  defect  shape.  By  minimizing  the  objective  function 
F  with  respective  to  the  parameters  by  any  of  the  optimization 
methods,  the  characteristics  of  the  defect  can  be  estimated. 

The  computational  process  in  inverse  problem  solution  is  shown  in 
Fig.10.  It  requires  solving  for  the  vector  of  design  parameters  h.  We 
first  generate  the  mesh  from  the  latest  parameter  set  h  to  get  the 
corresponding  finite  element  solution  for  A  and  then  compute  B. 


Fig.  8.  Voltage  induced  on  pickup  coil 
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From  B  we  evaluate  the  objective  function  F.  The  method  of 
optimization  used  will  dictate  how  the  parameter  set  of  device 
description  h  is  to  be  changed  depending  on  the  computed  F. 

Zeroth  order  optimization  is  practicable  in  terms  of  avoiding  the 
horrendous  programming  complexity  of  first  order  optimization 
although  computations  are  extensive  [14].  The  Genetic  Algorithm 
(GA)  which  is  a  zeroth  order  optimization  method  is  good  at  handling 
potentially  huge  search  spaces. 

Its  fitness  score  f  is  defined  in  terms  of  the  object  function  F. 
Although  our  object  function  F  as  defined  in  (5)  is  to  be  minimized, 
the  fitness  score  f  has  to  be  maximized  for  the  genetic  algorithm.  We 
therefore  define  the  fitness  score 


(6) 

According  to  the  methodology  of  optimization  using  GA  as  shown  in 
Fig.  11.  first  we  randomly  generate  hundreds  of  vectors  h  (each  called 
a  chromosome)  and  this  set  is  termed  the  initial  population.  With 
parallels  to  evolution,  a  new  generation  is  to  be  created  based  on  the 
best  of  this  population. 


Fig.  10.  Design  Cycle  for  the  computation  process 


Fig.  1 1 .  Optimization  Using  the  Genetic  Algorithm 

Then  the  fitness  score  for  each  h  is  calculated  and  checked  as  to 
whether  there  is  a  score  at  1  or  close  enough  for  our  purposes.  This 
computation  involves  computing  F  according  to  (5)  and  therefore  a 
finite  element  solution  for  that  h .  If  there  is  no  f  satisfactorily  close  to 
1 ,  then  the  design  parameters  are  changed  according  to  the  GA's 


classical  way  of  selection,  namely  selection,  crossover,  and  mutation 
Another  reason  for  using  GA  is  that  it  is  also  inherently  parallel  so  that 
it  may  be  easily  adapted  to  computations  on  the  graphics  processing 
unit  [15]  Though  GA  is  practicable  and  gives  a  better  solution,  it  is 
slow  when  compared  with  the  gradient  optimization  methods. 
Therefore  NVidia  GPU  parallel  computing  architecture  can  be  used  to 
solve  our  problem.  In  our  normal  programming  (single  CPU),  the 
fitness  value  is  calculated  for  each  chromosome  h  one  by  one.  When 
the  population  is  high  it  takes  a  very  long  time  to  converge.  Therefore 
we  launched  kernels  on  fitness  value  calculation.  For  that  we 
launched  GPU  threads  and  blocks  of  the  same  number  as  the 
population  size  (N)  so  that  the  fitness  value  will  be  calculated 
simultaneously  for  each  chromosome  in  the  population  as  shown  in 
Fig.12.  Therefore  in  this  paper,  we  evaluate  the  proposed  algorithm 
by  applying  the  GA  to  a  numerical  NDE  problem. 


Fig.12.  The  Parallelized  Process  of  the  GPU  [15] 


NUMERICAL  MODEL 


The  numerical  example  in  Fig.13.  is  used  to  validate  the  proposed 
algorithm.  The  coil  (with  jur  =  1.0,  and  current  density  J  =  ±5  x  104 
A/m2)  excites  the  magnetic  field  in  the  steel  plate  (with  /ur  =  100.0 
and  current  density  J  =  0.0).  The  conductor  is  surrounded  by  air 
(with  jur  =  1.0  and  current  density  J  =  0.0).  The  magnetic  field 
density  in  the  y  direction  By  is  measured  at  y  =  4.5  cm,  8  cm  <  x  < 
12  cm  using  10  points  in  the  interval  as  shown  in  Fig.13.  labeled  as 
the  measuring  line. 


A=0.0  J=  -50  kA/nv  Coil  J=50kA/m: 

Air 

Hr  =1.0 

J=0.0 

.  Measuring  line 


4 

A=0.0 


Steel  Defect 

Jlr =100  Z31(Air) 

J=0.0 


A 


X 


v 


<- . -> 

20  cm 

Fig.13.  Numerical  Model  for  defect  characterization 
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On  each  node  on  the  defect,  the  vertical  displacements  are  selected  as 
design  parameters.  In  our  numerical  model  we  have  8  geometric 

parameters  contained  in  the  vector  ^  (Y1*  V2'  V3*  V4’  V5»  V6»  V7’  V8). 
The  measuring  line  located  at  y  =  4.5  cm,  is  sampled  into  10  equally 
spaced  points  and  tolerance  boundaries  are  set  to  0.5  cm  <  h  <  3.5 
cm.  Each  design  variable  is  represented  by  10  bits.  For  testing  we 
took  a  particular  defect  h  as  h  =  {2  1  -2'  2-4,  2.2, 2.7,  3.0.  3.2,  3.5} 

cm  and  computed  the  field  (. BMeasured ). 
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Fig.  14.  Optimum  shape  of  the  Reconstructed  Defect 

Now  our  algorithm  has  to  reconstruct  h  to  match  the  “measurements” 
(B Measured)-  Design  parameters  are  changed  in  every  iteration  and  the 
infinite  element  mesh  newly  generated.  After  computing  A  by  finite 
elements  the  magnetic  field  (B Calculated)  is  computed  and  the  objective 
function  F  is  evaluated.  When  the  object  function  F  is  minimum  the 
iterations  are  stopped  and  h  is  found.  Fig.14.  shows  the  optimum 
shape  of  the  defect  after  200  iterations  for  a  population  size  of  200. 
This  shows  a  70%  accurate  reconstruction  of  the  defect.  This  error  is 
computed  from  the  equipotential  lines  for  the  finite  element  solution 
of  the  magnetic  vector  potential  Fig.15. 


Fig.  15.  Solutions  in  Equipotential  lines  for  the  Numerical  Model 

SMOOTH-SHAPED  DEFECT 

In  inverse  problem  design  optimization,  getting  a  practically 
manufacturable  shape  is  important.  An  erratic  undulating  shape  with 
sharp  edges  arose  when  Pironneau  optimized  a  pole  face  to  achieve  a 
constant  magnetic  flux  density  [16].  Their  results  are  shown  in 


FigJ_6.  The  nonsmooth  jagged  contour  in  FigJ_6b.  that  they  realized 
is  practically  not  a  manufacturable  shape.  This  they  addressed  by 
smoothening  the  pole  face  as  in  Fig,16c  manually. 


"AY  i />£^i§v 

llos 

_ ^  l/ffii 

a:  Start  b:  Optimized  c:  Smoothened 


Fig.  16.  Jagged  Pole  Face  of  Right  Half  of  Recording  Head  [17]. 

In  our  more  generalized  numerical  shape  synthesis,  the  problem  of 
jagged  shapes  was  overcome  by  imposing  constraints  [18].  In  a 
subsequent  paper  to  smoothen  a  surface  we  took  the  final  result  and 
set  the  coordinates  of  each  point  on  a  surface  being  shaped  to  the 
average  of  that  node's  coordinates  and  those  to  its  left  and  right  [191. 
We  also  used  element  by  element  matrix  solution  to  speed  up  the 
solution  process  [201. 

Since  our  design  optimization  is  for  defect  characterization,  there  is 
no  need  to  impose  constraints  to  get  a  smooth  manufacturable  shape. 
But  we  have  to  get  a  single  defect  with  a  realistic  shape  to  estimate 
whether  or  not  to  pull  the  vehicle  bearing  that  defective  hull  out  of 
service.  In  this  instance  however,  so  as  to  maintain  a  realistic  shape 
with  a  single  defect  we  imposed  the  constraints  as  hg>h1;  hy  >h2;  h6 
>h3;  and  h5  >h4;  Here  in  the  inequality  the  left  term  represents  the 
upper  surface  of  the  postulated  defect  and  the  right  term  the  lower 
surface.  In  effect  these  inequalities  ensure  that  the  surfaces  do  not 
cross  each  other.  Therefore  we  could  get  single  and  realistic  defect  as 
shown  in  Fig.15. 


CONCLUSION 

According  to  our  investigation  when  the  depth  of  the  defect  increases, 
it  is  hard  to  detect.  We  also  found  that  when  we  compare  horizontal 
and  vertical  defects,  that  a  horizontal  defect  has  a  higher  chance  of 
being  detected.  While  these  results  might  be  intuitive,  a  surprising 
finding  is  that  our  results  showed  that  when  a  defect  is  at  an  angle  of 
45°  from  the  surface,  it  has  a  higher  chance  of  being  detected. 

This  paper  also  presents  a  finite  element  technique  for  solving  inverse 
problems  in  magnetostatic  NDE.  Defect  shape  reconstructing  using 
the  genetic  algorithm  optimization  method  is  presented  and  validated 
using  a  numerical  model.  We  also  imposed  constraints  in  the  system 
to  get  a  realistic  single  defect  reconstruction  that  is  smooth. 
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